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Abstract 

A procedure to solve few-body problems which is based on an expansion over a small parameter 

is developed. The parameter is the ratio of potential energy to kinetic energy in the subspace 

. . , of states having not small hyperspherical quantum numbers, K > Kq. Dynamic equations are 
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subspace with K > Kq is taken into account in a closed form, i.e. without an expansion over basis 



. functions. 
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I. INTRODUCTION 



Below an approach to solving few-body problems which is based on an expansion over 
a small parameter is developed. The parameter is the ratio of potential energy to kinetic 
energy for the states with hyperspherical numbers K exceeding some limiting value Kn. 
Roughly speaking, the parameter is K^"^ . The method is a development of that of Ref. j^. 
An expansion over the parameter K^'^ has been given there for solving large systems of 
linear equations that arise in bound-state problems in the framework of the hyperspherical- 
hyperradial expansion, [ij] The method 2] is efficient for this purpose {4, 5]. However, 
for A>3 it is the calculating of matrix elements entering those systems of equations that 
requires a massive computational effort. The difficulty stems from a swift rise of a number 
of hyperspherical states with the same K a,s K increases, or a number of particles increases. 
Selection of hyperspherical states to reduce the effort, see , is efficient for A=3 and 

4 bound-state problems only. Such a selection is not justified in reaction calculations, in 
particular. The problem is removed in the method below since no expansion over basis states 
is employed here for K > Kq. 

Recently a considerable progress in methods for solving few-body problems has been 
achieved. However, those developments have limitations, and the latter are removed in 
the present method. In particular, the well-known Green Function Monte Carlo method 
to be mentioned in this connection is the method to calculate a bound state of a system, 
and it is not suit to calculate reactions. (Although the simplest scattering problems may 
be considered in it frames.) Besides, this method is not convenient in the respect that it 
provides separate observables, such as an energy or a size, as a result of a calculation but it 
does not provide the wave function of a bound state that could be employed in subsequent 
calculations. Unlike this method, the method below is suitable for calculating reactions of 
a general type. And when in its frames one needs to use a bound state wave function one 
need not recalculate it completely each time. 

Recently a way was found to extend the Faddeev-Yakubovsky A=4 calculations over 
the energy range above the four-body breakup threshold [8]. However, Yakubovsky type 
calculations require too much numerical effort even in the A=4 case. Amount of calculations 
is considerably less in the scheme below. 

At solving few-body problems with expansion methods convergence of expansions for 
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calculated quantities was accelerated with the help of the effective interaction approaches. 
Such approaches were developed in the framework of the oscillator expansion [9] and the 
hyperspherical expansion [l^ . In their framework a true Hamiltonian is replaced with some 
effective Hamiltonian acting in a subspace of only low excitations. When, formally, the latter 
subspace is enlarged up to coincidence with the total space an effective Hamiltonian turns to 
a true one. An effective Hamiltonian is constructed from a requirement that its ingredients, 
as defined in a subspace of low excitations, reproduce some properties of the corresponding 
ingredients of a true Hamiltonian in the total space. It has been shown jo], [l^ that this, 
indeed, leads to an improvement of convergence of observables considered. 

Higher excitations are disregarded in such type calculations. It is clear, however, that 
correlation effects related to higher excitations cannot be reproduced by any state vector 
lying in an allowed subspace of only low excitations. For example, let us consider the mean 
value, (^^ol-f^l^o)) of such an "observable" as a true Hamiltonian. It follows from the varia- 
tional principle that an approximate state \E'o supplied with such a method provides poorer 
approximation to the true (\I'o|-f^|\l'o) value than \I'o obtained by the simple diagonalization 
of a Hamiltonian in the same subspace of low excitations. And even the value of (\E'o|-ff|^'o) 
obtained with the latter \E'o is a very poor approximation for realistic Hamiltonians. On the 
contrary, the method given below provides an approximate state vector that is apparently 
close to a true state vector both as to its low excitation component and its high excitation 
component. 

And speaking of reaction calculations in the framework of Eq. (1401) below, [H — a)'^ = q, 
one should in addition take into account that a rate of convergence is determined not only by 
properties of the Hamiltonian H but also by those of the source-term q. But these properties 
are apparently ignored at constructing effective Hamiltonians. Therefore one cannot expect 
fast convergence in all the cases, especially for source-terms q corresponding to strong- 
interaction induced reactions. On the contrary, the method described below provides state 
vectors genuinely close to the true ones both for bound state problems and any reaction 
problems. 

In the next section the bound state case is considered. In Sec. 3 modifications to treat 
reactions are listed and a numerical estimate of the rate of convergence of the method is 
done. Some comments on computational aspects contain in Sec. 4. 
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II. BOUND STATES 



We consider the eigenvalue problem 

{H-E^)^^ = 0, (1) 

where H = T + V is an A-body Hamiltonian. We split the whole space of states into 
the subspaces with K < Kq and K > Kq and we denote and the components of 
the solution that lie, respectively, in these subspaces. At a proper choice of Kq kinetic 
energy T of a state belonging to the second of the subspaces is much larger than its potential 
energy. Indeed, 



T = Tp + - ^, 

where is the hyperangular momentum operator acting on a hypersphere, p is the hyper- 



radius, and Tp is the hyperradial energy operator. The eigenvalues of the K operator are 
K{K + n — 2) where n = 3A-3 is the dimension of a problem. Thus (T) is large for states 
having large K and not too large space extension. We choose Kq in a way that for K > Kq 
one has, in a rough sense, 

K^ 



»|[(V + T,-E,)vl/,]^|. (2) 



Here [. . ]k denotes a component of a state with a given K. Eq. ([2]) is to be fulfilled for all 
configurations that contribute significantly to a solution. The corresponding p values range 
within the configuration space extension of a solution. These p values are such that p^ is 
less than, or about, A{r'^), where r is the single-particle size of a system. 

At these conditions one may express the component ^\ in terms of ^\ perturbatively and 
obtain equations for the latter component alone. Let us define projectors onto the K < Kq 
and K > Kq subspaces as Pkq and Qko, respectively. Let us express formally in terms 
of 

ml = ~TK,{Ey:)Vm\, (3) 

where 

VK,{E) = [QK,{H-E)QK,r' (4) 

is the Green function defined in the second subspace. It is taken into account in ([3]) that 
kinetic energy is diagonal with respect to K. It is convenient to define Vkq as acting in the 
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whole space and to rewrite it in the form 



Qko [Qko{H — E)Qko] ^ Qko- 



(5) 



Substituting Eq. into the relationship 



(6) 



one gets the equation for \E'^ alone, 



(7) 



The quantity PkqVT KoiEx)V represents the exact effective interaction arising due to 
coupling of the complementary K > Kq subspace to the K < Kq subspace. 

We shall solve Eq. (I7j) perturbatively. We write in ^ H — E = L + U and we use an 
expansion 



'^Ko{E) = Gkq — GkoU{E)Gko + G kJJ {E)G {E)G Kq — 



(8) 



where Gko = Q KoiQ k^LQ Kq) ^Qkq- With the choices of L below it has no non-zero matrix 
elements between the subspaces with K < Kq and K > Kq, i.e. [Qkq, L] = 0. Then 

Gko = QkqL ^ = L ^Qkq- 



For performing calculations in the coordinate representation we choose Gko as follows, 

-1 



G 



2m p 



— +Wip) 



Qko- 



(9) 



It is convenient to represent ([9]) as a sum of contributions from various K values. 



K>Ko 



(10) 



Then 



K{K + n-2) 



2m 



+ W{p) 



-1 



S{p - p') 



p 



:Jr^Y.^KumKAO- (11) 



Here ^ and ^' are n-dimensional space vectors, W{p) is a subsidiary interaction, ^ = ^/p, 
^' = C,'/ p, and Yku form a complete set of orthonormalized hyperspherical harmonics having 
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the same K. The hyperangular factor entering here may be represented with the simple 



expression (e.g. 



ul) 



T.'^LmK.ii') = ^^r(7)cx(e ■ (12) 

V 

where C^(x) is the Gegenbauer polynomial. 

The choice ( fTTll of qk is done to facilitate Monte-Carlo calculations of matrix elements. 
At this choice one has in (JHl) 

^i\v{E)\o = m\o + 5(e - a [T,-E- wip)] , (13) 

T.^-^(^^"-^4:)- (14) 



2m \(ip^ p dp 

To perform calculations in the momentum representation we suggest the expansion 
with a modified Gko, Gko = [nV(2m) - Eq + W(J1)]-^Qko, 



{tt\Gk^\7i') 



nV(2m) -Eo + W(U) 
6{n - ttQ - n-("-^)^(n - nO(2 ■ n-/')-^ Y.k<kSK + 7)C^(vr ■ ttQ 

nV(2m) - + W{Ii) 

_ n-("-i)5(n - n')(2 ■ n-/')-^ Ek>k^k + i)ci{n . V) 



(15) 



nv(2m) - ^0 + w{n) 

{tt\U{E)\tt') = (7r|\/|7r') - 6(71 - 7r')[E - Eq + W(n)]. (16) 

Here vf and vf' are n-dimensional momentum vectors, 11 = |7f|, tt = vf/II, tt' = vf'/II, and 
Vr(n) is a subsidiary interaction. The quantity Eq is a fixed energy chosen to be close to 
E\ sought for. 

Roughly speaking, the expansion goes over K^"^. As Kq increases relative contributions 
to a solution from subsequent terms in the expansion ([8]) decrease. Taking Kq sufficiently 
large we retain only the lower terms in the expansion. 

The subsidiary interaction W{p) ~ V{p) or H^(n) ^ ^^(n) is intended to accelerate 
convergence of observables of interest when Kq increases. A better choice of subsidiary 
interactions would be such that they include spin-isospin operators. Let us suppose that 
calculations are performed in the coordinate representation. For a conventional NN inter- 
action that includes static local central and tensor components Vioc plus components that 
depend on angular and linear momentum a possible good choice is the following. Let us 
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consider the K = component in the expansion of Vioc over hyperspherical harmonics. This 
component is the resuh of averaging Vioc over a hypersphere. It has the structure F{p)0 
where O = "^Oiij) is an operator that depends on spin-isospin variables. The operator O 
is symmetric with respect to particle permutations. Therefore it may be represented as 

where / labels irreducible representations of the permutation group of A particles, /i labels 
basis vectors belonging to a representation [/], {9jl } is the corresponding orthonormalized 
set of basis functions, and 0/ is defined as follows, 

{9^^^\d\9^;)') = Sff,6,,,Of. (17) 

We then choose W as 

w = F{p)J20fJ2\elf^){e^^^\. 

f M 

At this choice, V and W cancel each other to a large degree in the difference V — W entering 
U. This allows employing a smaller Kq value. The Green function Gkq = X^i^ 9k becomes 

«l9A-l«'> = ^^^|^^rwcj({.r) 



X 



The quantities Of may also be varied around their values from ( |T71) . To simplify the pre- 
sentation we did not include a spin-isospin dependence in the formulas above. 
We set in ([7]) 

*i = E^?"^' E, = J2Et\ (19) 

n n 

where v|/^(") and E^^^ correspond to the n-th order in the expansion over GkqU in ([H]). We 
then get from ([7D, dHD 

Pxo(i^-i5^f)*r = 0, (20) 
Pk,{H - Ef^)^T = + PK,VGK,V^>f\ (21) 

+PK,VGK,V^f\ (22) 



If ^ is a solution to Eq. fl^Tj) then \1>^^^^ + c\I/^^°^ with an arbitrary c is also a solution. 

The same holds true as to ^^[^^ in (I22D- To get a unique solution it is sufficient to impose 
the normalization condition 

(^ii*i) = (^ri^r)- (23) 

This gives in the ffist and second order, respectively, 

(vl/f |vl;f)) + (vl/f)|vl/f)) = 0, (24) 

i^f) + (*ri*f ) + (^ri^r) = o. (25) 



Taking into account time reversal invariance of the operators entering ( ETj) . (1221) it is seen 
that the matrix elements in (1231) and (1251) are real. Therefore (1211) and (l25l) turn to 



(Vf)|^f)) = 0, (26) 

(Vf )|^f + (l/2)(^f = 0. (27) 

Taking scalar products of Eq. (1211) and Eq. (12^ with and making use of Eq. 
we obtain, respectively. 



77.1.1; _ \ A I ' ^ -t^O ' I ^ X I /r,Q\ 



To get (j29l) Eq. (l26l) was employed. 

We seek for the component as an expansion over the hyperspherical basis. In the 
coordinate representation, 

(30) 

K<Ko;i' 

Here a^i and r^j are particle spin-isospin variables, F^u are basis functions that we consider 
to be orthonormalized. They are combinations of basis hyperspherical harmonics and basis 
spin-isospin functions. It is implied here and below that all the summations over K include 
only K values of a given parity. Let us write down similar expansions for \[''*^"^ 



A ' 



K<Ko;u 
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so that xku = ^nXxl- -^1^- (SO]), turn into equations for the expansion coefficients 

(n) 

(^x-i?f)x£+ E (^HV^I^V)xS?l, = 0. (31) 

K'<Koy 

(TK-Ef\^il+ (^HV^I^V)xiJl, = (irHV^G^o^l^r) + 4'^xi?l. (32) 

K'<Ko-y 

K'<Ko;u' 

= {Ku\VGK,V\^>f^) - {Ku\VGK,,U{Ef^)GK,y\^>f^) + + Ef\^^^. (33) 

Here denotes the hyperradial operator of kinetic energy, 

K{K + n-2 
2m p 

In the notation above {Ku] . . .) = {Fku\ ■ ■ ■) and {KulVlK'u') = {Eku\V\Fk'u')- These 
quantities are defined in a obvious way. We recall that the equations written down include 
K values only within a finite range, K < Kq. The zero order equations fl3ip are the 
standard ones that arise when coupling to states with K > Kq is disregarded. The higher 
order equations just take this coupling into account. 
Eq. ([21] reads as 

[p--'dp Y xPMxil{p) = 0. (35) 
The condition ( l35l) is to be added to Eqs. ( l32l) . Let us suppose that Eqs. ( 13T|) and (132|1 



K=T, + —^—, (34) 



are solved via an expansion of x^ku ^"^^ Xku o'^^^ same hyperradial basis with the same 
number of basis functions retained. The linear equations arising in this case from Eqs. (132|) 
are linearly dependent. In general, one should remove one of these equations and replace it 
with the hnear equation to which Eq. turns. Eq. (^3]) becomes 

I P^'-'dp J2 xtip)x^iiip) + im[x^il{p)]' = 0. (36) 

K'<Ko;u' 

This should be used similar to Eq. ( 135|) . If instead of ( l30l) a hyperspherical expansion is 
employed within a momentum representation calculation similar equations may be written 
down proceeding from fl2ni) - fl22l) . 

The complementary K > Kq component \I'^ of a state sought for may be written as 

= E ^^^"^ (37) 



where vl/^*^"-' signifies a contribution having the n-th order in GkqU, and v[/^'-°'' = 0. Then 
one has 

v^f ^ = -GKoV^f\ (38) 

) = -G^oV^^f ^ + GKoU{E^^^)GKoV^f\ (39) 

The component has been obtained above in the form of a hyperspherical expansion. 
Therefore one may store it and use in various apphcations. The complementary component 
then may be reconstructed as a simple quadrature fl38|) . fl39l) . 

If, for example, \1/a is calculated up to the n > 1 corrections, i.e. vEf"^^^ = vl/^^*^) -|- v]/^'^^) -|- 
then the average energy = (^I^PP^^I^f differs from the exact 



value in terms only of the third order and higher in the expansion over Gkq- In particular, 
the second order energy ( l29l) is correctly reproduced with Ex- Indeed, according to the 
variational principle the difference between Ex and the exact Ex value includes the term 
(^^aI^^'I^^a) and powers of the term (5^a|5^a)- Here 5^a = ^f"'* - ^f^''- We have 
5^ A ~ {GkoY while presence of H in the above matrix element changes the net power in 
Gko from (Gko)^ to {GkoY- 

Basing on Table 4 in Ref. ^ one infers the following. When only the above considered 
n = 1 correction is retained the choice Kq = 14 ensures the correct binding energy at the 
accuracy level better than 0.1 MeV in the A=4 bound state problem with a realistic NN 
interaction that includes a strong core. The net number of HH with i^' < 14 entering the 
problem does not exceeds several hundreds which is acceptable. 15(] 



III. REACTIONS 



1 . We consider a dynamic equation of the form 

{H - a)-^ = q. (40) 

Here a is a subsidiary complex energy, and g is a given state. Reaction amplitudes may 
be obtained from ^(o") in a simple way, see e.g. the review \x^- The approach extensively 
applied to perturbation induced reactions and proved to be very efficient. Any strong- 
interaction induced reactions can also be treated with this approach. 

The solution \E' is localized. Therefore the procedure quite similar to that described above 
is applicable also here. One represents as + and obtains these components as sums 
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of successive approximations, = \|/'(")^ v]/^ = -qjKn)^ where the meaning of notation 
is the same as above. One has 

Pk,{H - a)¥^'>'^ = PK,q, (41) 

Pk,{H - a)^'(i) = Pko{VGk,v¥^'^ - VCkoQ), (42) 

Pk,{H - (t)^'(2) = Pk, [-VGK,U{a)GKoV^'^'^ + VGk,v¥^'^ + VGK,U{a)GKoq] • (43) 

As above these equations may be rewritten as coupled equations for coefficients of the hy- 
perspherical expansion in the coordinate or momentum representation. The complementary 
components are obtained from as quadratures, 

^^^^^ = -GKoV^'^''^ + GKoq, (44) 

= -Gk,V¥^'^ + GKoUia)GK,V^'^'^ - GK,U{a)GK,q. (45) 

When, for example, it is sufficient to account for only the n = 1 corrections in $(0") = 
(^(cr)|^((T)) one need not calculate the components. 

To estimate roughly the required Kq value we note that the large-distance decay of \P 
in the configuration space is determined by the imaginary part of the wave vector k = 
[{2m/h^YaY/'^ . Let us write a = aR + ai and denote R = (ImA;)^-*^. Let us suppose that 
a calculation is performed in the coordinate representation, and the expressions ( JTOl) . (|TT1) . 
([12]) for Gko and (|T^ for U with E = a are used. Then similar to ^ one may estimate the 
required Kq value from the condition 



2,nflH''° + — J (46) 

A typical o"/ value is 10 MeV, and a required range of ctr values is about the same as a 
range of energies considered in a problem. When aR is not too high Eq. fj46l) is fulfilled for 
acceptably low Kq values. (We shall not discuss the point on a precise V value to be put 
there.) 

When, however, the quantity gr is high the expansion ([8]) of the Green function converges 
quickly only for large Kq values. (The deceleration of convergence is caused by both terms 
Tp and a in f/, while they may compensate each other only in part.) To speed up the 
convergence, one could remove the contribution Tp — a from U and to account for it in 
qk- However, this would hamper Monte-Carlo integrations because of the rapidly changing 
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hyperradial Bessel and Hankel functions that would enter in this case. One may avoid 
these comphcations if one performs calculations in the momentum representation. In this 
case one uses Eq. (|T5l) for Gko with Eq = a, 



{n\GKo\Tr ) 



_ _ 6{n - n') - n-("-i)5(n - n')(2 • n-/')-^ EK<KoiK + ^)Clin ■ n') 



n7(2m) -(T + vr(n) 

„ (K ^-v]C7Jrr ■ n') 

(47) 



n-("-i)5(n - n')(2 ■ vr«/2)-i J:^^^^{K + j)Cl{n . n') 



nV(2m) - a + w(n) 
Correspondingly, in Eq. ( |T6l) for U one replaces E — Eq with zero, 

{7r\U{E)\n) = {7i\V\n) - 6{7v - n)W(n). (48) 

In this case the condition 

Kq + ^— > V 



is sufficient for quick convergence of the expansion ([8]) for the Green function. Considering 
the role of subsequent terms in ([8]) in this case one should take into account that if a 
coordinate representation wave function f{^) is localized within a hyperradius p then the 
momentum representation quantities (y/^j,(7r)|/(7r)) are very small at the 11 values such that 
Up <^ K + {n — 2)/2. (Irrespective to the mentioned condition, the condition cr^ >> V 
also leads to quick convergence of the expansion ([H]).) While Eqs. ( H71) . (HSl) are required 
for performing calculations that involve high values, these equations, of course, may be 
employed at low cr/j values as well. 

Let us perform a rough estimate of efficiency of the latter version of the approach. Let 
us consider the A=4 case and adopt the Kq value equal 14. Let us estimate the relative 
role of the correction '^^'^\a) with respect to ^^^\a). For this purpose let us compare the 
contributions of these corrections to the net transform (see jl2'|) ^{cr) that correspond to 
K = 16. These contributions are ^^'ilwi^) = (^iJLiel^iJlie) and $iJL^i?^(a) = {^kIw + 
^K=w\^K=m + ^K=i6)- We take in (HTD W{Il) = and perform the calculation in the 
coordinate representation. For estimate purposes we can assume that \E'' is given, and with 
its help and ^'^(^^ are subsequently calculated as 

^''^'^ = -GK,V¥ + GKoq, 

^h(2) ^ (49) 

(These expressions are not the same as fl44p . since we consider \E'' to be known here.) Let 
X^Kuip) and x^Kuip) tie the coefficients of expansions of i^'*^^) and i^^(^) over hyperspherical 
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harmonics. Then Eq. (HUj) turns to 



JO 

where 

and the free motion Green function is 

gPip^p') = {2m/n'){i7,/2){ppr'JK+,{<yp<)H^K\,{<yp>)- 

Here p< = min(p, p'), p> = max(p, p'), 7 = (ra — 2)/2. One may also write at the K = Ko + 2 
value 

{V^''^'^)kAp')= E Wv(p')xiJl(p')^ W.(p')xS(p')^^o,o(p')xS(p')- (50) 

In what follows we omit the subscript u and perform the estimate up to multiplicities in u 
both in {^ilj^ilj and le + ^^kUMU + ^kIJ- Thus we use 

POD 

X^kIM= / 9P=ieiP,P)Vo,o{pViU{p')iprdp'. 
Jo 

For Xk^=w{p) we adopt the model 

1 pikp 



where ctq = 50 MeV, and {TikY /{2m) = a = + icr/.[16[ The expression (ISTjl ensures the 
correct asymptotics at large p values. We set cr/ = 10 MeV that is a good choice to invert 
the transform, and po = 2 fm. We take V = ^(^ij) we employ the Gaussian potential 
V{r) = Vo exp(— with the parameters Vq =67 MeV and b =1.5 fm that corresponds 
to a triplet potential reproducing the scattering length and the effective range. The values 
of $/<Li6(^) $^'|j^g^'' (a) are shown in Fig. 1 as functions of a^. It is seen that the 
second order correction to is, indeed, of minor importance as compared to the first order 
correction. 

2. Let us also consider reactions at low energy when only two-fragment channels are 
open. Suppose that the dynamic equation {H — E)^i = for continuum spectrum states 
"^i is solved and channels are open. One may use the well-known ansatz 

N 

i=i 
13 
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FIG. 1: The relative effect of the first and the second order corrections to the Lorentz transform 
$((t) = (^'1^'). One has a = cr/j + id/, and the o"/ value is taken to be 10 MeV. 



where (f)i^ and (f)f^ represent the "channel" states of two possible types, while X is localized 
and is sought for as an expansion over hyperspherical harmonics. Let K values up to some 
Kmax are retained in the expansion, and PRmax is the projection operator onto the subspace 
of those harmonics. The corresponding approximate equations include those to determine 
the expansion coefficients Xku{p) at reaction amplitudes fij supposed to be "given". These 
equations may be written as 

Pk^.AH - E)X = PK^^^q, (53) 

where 

N 

i=i 



To fix fij one adds linear equations. 



li 



We note that Eqs. ( l53l) may be solved perturbatively similar to Eq. ( HOl) . If < is 
not too close to the three-fragment reaction threshold so that p = {X\p\X) is not too large 
and the inequality (H6l) is fulfilled for moderate Kq values then one can employ Kq <^ Kmax 
and retain only lower terms in the corresponding expansion over GkqU. (Of course, this 
procedure is of limited use in the case of n-d scattering.) 
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IV. COMMENTS 



1. If the partial waves Xku{p) from (IHUj) are sought for as an expansion over a set of 
hyperradial functions the above dynamic equations turn to systems of hnear equations. 
These systems may be of a large size. Then they can efficiently be solved with a version of 
the method of Ref. (1^, i.e. again using an expansion over another parameter of the Kq^ 
type. We note that when this method is applied subsequent iterations are identical to each 
other in their form so that it is easy to perform a required number of them. (Their number 
may typically be about ten or so to provide an accurate solution Q].) On the contrary, in 
the method described here an increase in a number of iterations means an increase in the 
dimension of the corresponding integrals. Therefore the present method is practical only 
when low order corrections are sufficient. This is the price for the elimination when K > Kq 
of an expansion over basis states in the present method. 

2. In the case of reaction calculations one passes from solutions of Eqs. fl4Up to reaction 
amplitudes as follows (see e.g. [3]). The quantities of the $(cr) = {^{cr)\^{a)) type are 
formed. These quantities are integral transforms of response-like form factors that determine 
reaction amplitudes. So, to pass to reaction amplitudes these integral transforms are to be 
inverted. To perform a satisfactory inversion of the transform $(o") one needs to use its 
values in rather many cr points. But one need not solve Eqs. fHUj) for the corresponding 
many a values employed. A better approach is to solve these equations for rather a small 
number of a values and to obtain $(ct) for a larger set of cr values via interpolation. The 
transforms $(cr) are smooth functions and this procedure is safe and accurate. 

3. Those matrix elements in the above equations which are related to contributions from 
K > Kq are to be calculated with the Monte-Carlo method. It still should be verified 
numerically whether the Monte-Carlo integration is efficient enough for this purpose. The 
existing experience testifies to that the Monte-Carlo integration is suitable at least in the 
case of matrix elements entering Eqs. (I3T1) even when K values are rather large (see e.g. 

The coordinate representation matrix elements above that correspond to the n = 1 
correction have the structure 



K>Ko 



K(K ^r, —0\ 1^ 



2m 

(54) 
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Here, for example, F2{pi') = -^2(^0 = ('C'I^^a'''')- The notation or d^' refers to hyper- 
angular integrations. One deals with similar type integrals also in calculations that involve 
Eq. flHHl) . Apart from a direct Monte-Carlo integration, in some cases it is expedient to 
take the argument of the Gegenbauer polynomial in (15^ as a new variable, see Appendix, 
and to integrate over this variable with use of the regular Gauss-Gegenbauer quadratures. 
While all other integrations are to be done with the Monte-Carlo method. This can also be 
done in the case of a momentum representation calculation. Last lines in Eq. (fT^ or (H7|) 
are to be used in this case. 

It is convenient to use permutational symmetry of states to simplify calculations of the 
n = 1 correction. For example, when one retains only a two-body force, V = X]j<j ^(^j)) 
one can write 

A{A- 



{^i\VgV\^2) 



{^i\V{12)gV{12)\^2) + 2{A - 2){<i!i\V{12)gV{13)\^2) 

■ (55) 



+ ^\ ^,\V{12)gV{34)\^il2) 



When a three-body force is retained similar relationships could be written as well. Eq. ([51 
is written up to spin-isospin variables. When Green functions Gkq or g from flSSp are spin- 
independent it is convenient to include the intermediate spin-isospin factor = I 
in them, where {6^} is a complete set of spin-isospin states, c.f. (|TS|) . 

4. Let us comment on the n = 2 correction. Suppose that a conventional NN interaction is 
employed that includes local central and tensor components plus components depending on 
orbital and linear momentum. Contributions from local components of such an interaction to 
the n = 2 correction have the following structure in the case of the coordinate representation 
calculation, 

^ ^ (5 + 7)(i^' + 7)r^(7) 



(2 ■ W2)2 
K>Ko K'>Ko ^ ' 



X / p--Hpdi^di2diz 



r K{K + n - 2) 
2m 



+ W{p) 



K'{K' + n-2) 
2m p^ 



+ W{p) 



When n = 2 corrections are retained in a calculation sufficient accuracy is provided already 
with rather small Kq values. Then contributions to fl3Ul) only from not large K and K' 
are significant which facilitates the Monte-Carlo integration. It may also be noted that 
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contributions of (I56p type with the above mentioned non-local components of NN interaction 
include operators acting on the Gegenbauer polynomials. To disregard these contributions 
is a good approximation in many cases. 

5. As one could infer from Table 4 in [J] a Kq value required to ensure convergence is 
considerably smaller in the case of NN interaction with a super soft core than that in the 
case of NN interaction with a strong core. Therefore one probably could reduce a required 
Kq value also via transformation of dynamic equations to a form that involves a t matrix 
rather than an NN potential. (In the A=3 case such equations are the Faddeev integral 
equations but at A>3 there is no need to pass to the Yakubovsky type equations for this 
purpose.) 



APPENDIX 

When one takes ^ ■ ^' in (l54l) as a new integration variable one needs to define the whole 
set of integration variables in a way that the integrand remains smooth. This can be done 
e.g. as follows. Let us express ^ = {^i, . . . , in} in terms of another unit vector 57, 

n 

where gij is an orthogonal matrix such that its first column is gn = and gij is arbitrary 
otherwise. One then has i ■ i' = j 9ii9ijVj = Vi- Let us parametrize the components of fj 
as follows, 

fji = cos LP, fjj = Vj_i sin v?, j = 2, . . . , n, 

where Vi are components of a unit vector {) on a hypersphere in a n — 1-dimensional subspace. 
Taking into account that 

di = dr) = (sin Lp)'^~'^dvdLp 
one then may rewrite the integral ( |5^ as 



K>Ko 



E ^^r(^) / P'^-'dpdi'dvism^f^dLpF.ipi) 
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where the components of the n-dimensional unit vector ^ entering Fi are parametrized as 
follows, 

The integrations over dp, d^', and dv may be performed with the Monte-Carlo method while 
the remaining integration over d(p may be done with the help of regular quadratures. 
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